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Abstract. In this work, we introduce matrix- valued diffusion processes which de- 
scribe the non-equilibrium situation of the matrix models for the beta-Hermite and 
the beta-Laguerre ensembles. We also study the corresponding spectral measure 
process and empirical eigenvalue/singular value process with regard to their limit 
laws. 
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1. Introduction. 

The beta-ensembles have the physical interpretation as ensembles of one dimen- 
sional Coulomb gas in a neutralizing background with appropriate charge density 
[Dy,F], where /3 is the inverse temperature. For /3 = 1,2,4, matrix models of the 
beta-ensembles existed for a long time and are known as the *-orthogonal, *-unitary, 
and *-symplectic ensemble respectively, where * E {Gaussian, Laguerre, Jacobi}. 
(See for example, [Meh] and [Muir] and the references therein.) In recent years, as 
a result of the work of Dumitriu and Edelman [DE] , and that of Killip and Nenciu 
[KN], matrix models for the beta-ensembles with entries/parameters from classical 
distributions are now known for all values of /? > 0. 

In this work, we will restrict ourselves to considerations which are related to the 
matrix models of the /3-IIermite and the /3-Laguerre (Wishart) ensembles. In order 
to explain what we want to do, let us consider the matrix model of the /3-Hermite 
ensembles [DE], defined schematically by 
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where the entries on the diagonal and the subdiagonal of the n x n matrix J^ 
independent of one another. Prom the definition in (1.1), the joint density of the 
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independent entries of Jg is given by 



Wh{a, h) = Cnp JJ l>f_f} exp 
fc=i 




+ 2E^? 



i=l 



(1.2) 



where a^ = {Jp)ii, hi = (J^)i^i+i, and c„^ is a normalization constant whose value 
is given in (3.4). In this work, one of our motivating questions is the following: 
Is there a natural Markov process on Jacobi matrices for which Wh{ci-,h) is the 
stationary steady state density ? Our approach to this question can be described 
as follows. Since the entries on the diagonal and the subdiagonal of Jp in (1.1) 
are independent random variables, therefore, it is natural to seek a Markov process 
which is a product of statistically independent Markov processes on the diagonal 
and the subdiagonal of the matrix. In this way, the problem reduces to that of 
constructing one-dimensional Markov processes with the required properties. As 
it turns out, for the class of probability density functions p{x) satisfying Pearson's 
equation [P] (with support on an interval / C M with endpoints xi and X2, say), 
the problem of constructing a class of stationary Markov processes with transition 
probability density p{t,xo,x) for which 



lim p{t, xo,x) 



p{xo)pit, xo, x)dxo = p{x) 



(1.3) 



Xi 



has been the subject of investigation in [W]. The Pearson's family of probability 
density functions is quite broad, indeed it includes many of the continuous dis- 
tributions which we commonly use such as normal, Chi-square, beta, etc. For 
us, although the Chi distribution is not on the list, however, the solution of the 
problem for Chi-square suffices. This is because the Markov process correspond- 
ing to the pdf of the Chi-square distribution can be identified to be the square 
of a generalized Bessel process. Hence the correponding process for the Chi dis- 
tribution must be the generalized Bessel process itself and indeed, the transition 
probability density function for this process also satisfy the requirement in (1.3) 
above. (See Section 2.1 below.) As for the normal distribution with mean 0, it 
is quite well-known that the Ornstein-Uhlenbeck (OU) process describes the non- 
equilibrium situation [N]. In fact, one can readily understand Dyson's introduction 
of the matrix-valued Ornstein-Uhlenbeck processes in [Dy] from the perspective 
which we discussed above. Once the matrix-valued diffusion processes are in place, 
there are of course many questions one can ask about the statistics of the corre- 
sponding eigenvalue processes. In this work, we will begin with basic results on 
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the spectral measure and empirical eigenvalue processes. (In the case of the /3- 
Laguerre process, we will instead deal with the empirical singular value process.) 
Thus in a sense, we are trying to understand the non-equilibrium properties of the 
Coulomb gas for general values of /3 > and the these processes are basic. In this 
connection, we remark that for /? = 1,2,4, the eigenvalue processes induced from 
the corresponding /3-Hermite diffusion processes are actually different from those 
in [Dy], although there appears to be some subtle connections. We hope to return 
to this and other aspects of these processes in subsequent publications. 

The paper is organized as follows. In Section 2, for the convenience of the reader, 
we begin with a review of the generalized Bessel and other related processes that 
play a key role in this work. In Section 3, after a short description on the product of 
independent Markov processes, we introduce the /3-Hermite processes (J/3(t))t>o for 
general values of /3 > in the first subsection. Then we calculate the distribution of 
the eigenvalues and the first components of the normalized eigenvectors of Jp (t) for 
each t > 0. The latter quantity, as it turns out, is independent of t. For the spectral 
measure Ht = X^?=i /^j(^)'^A (t) associated with Jf3{t), we also give the distribution 
of ^ -^j^ f^jit) foi' each 1 < A; < n. In the second subsection, we consider the scaled 
process (Ji" (t))t>o, where JJ^ {t) = -^j^- The goal is to investigate the weak 
convergence (as n — t- oo), in probability, of the corresponding spectral measure fi^ 
and the empirical eigenvalue distribution i/" for each t > 0. Because of the well- 
known one-to-one correspondence between Jacobi matrices and spectral measures 
(see, for example, [Dl]), it is natural to consider /i". In this connection, we make 
use of the method of moments to show that jjl^ converges weakly, in probability, to 
a deterministic measure jit whose density function has the shape of a semicircle. In 
this way, we obtain a time-dependent semicircle law. Because the distribution of 
the first components of the normalized eigenvectors of Jp (t) is independent of t and 
the distribution of Y^j^i ^^ji't) is in fact identical to the one in [BNR], we can make 
use of the proof of Theorem 5.4 in [BNR] to conclude that d^p^n^, u^) converges in 
probability to as n tends to infinity. (Here dLP is the Levy-Prohorov metric.) As 
the authors are studying a rather different problem in [BNR] , it is quite remarkable 
that we can make use of their analysis here. In Section 4, we do the same for 
the /3-Laguerre processes and the associated /3-Wishart processes. In the case of 
the scaled /3-Wishart processes, the limiting law of both the spectral measure and 
the empirical eigenvalue distribution is that of a time-dependent Marchenko-Pastur 
law with a hard edge at and a soft upper edge at 4/9(t), where p{t) = 1 — e~*. 
Finally,for the /3-Laguerre processes, we show that the limit law for the empirical 
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singular value distribution is that of a time-dependent quarter circle law. 

2. The generalized Bessel processes. 

Let (b(t))t>o be standard Brownian motion on M, and 6 > 0. The Bessel process 
{R {t))t>o of dimension 6 starting from [RY] is the square root of the process 
defined by the SDE 

duit) = 2^/[[i^\db{t) + 5dt, u{0) = 0, (2.1) 

where the point is an instantaneous reflecting barrier for < 6 < 2 and is polar 
otherwise. When 5 is a positive integer, it is well-known that R^{t) is the modulus 
of a standard Brownian motion in M . For our purpose, recall that the the transition 
probability density function of R {t) is given by 

pf(xo,x)=t-i(^y 'xexp(-^i±^)/._,(^) forxo>0, (2.2) 



and 



pf(0,x)=2-Hit-lr('Q" .^-iexp(-^~), (2.3) 



where r(x) is the Gamma function, and /^ is the modified Bessel function of order 
V. Since we have i? (0) = with our definition, it follows that the pdf p^(x) of R (t) 
coincides with pf{0,x). We will also deal with the Ornstein-Uhlenbeck processes in 
this work. By definition, the one dimensional Ornstein-Uhlenbeck process with 
parameters (a, a) starting from is the solution of the SDE 

dv{t) = -av{t) dt + a db{t), v{0) = 0, (2.4) 

where a > and o" G M \ {0}. By the Dambis-Dubins-Schwarz theorem [RY], v{t) 
is a time-changed Brownian motion: 

/ 2at _ 1 \ 

y(t) = ae-^'w (^^^J > (2-5) 

where w{t) is a standard Brownian motion in M. Consequently, the random variables 
v{t) are Gaussian with mean zero and variance 

pW = ^(l-e-'"*), (2.6) 
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and the transition probability density function is given by 

Note that if v(t) = (fi(t), • • • , Vd(t)), where t'i(t), • • • , Vdif) are independent one 
dimensional Ornstein-Uhlenbeck processes starting from with the same parame- 
ters (a, o"), then from (2.5) above, we have 

Vi{t) = ae-'^'wi { ^^ j , i = 1, • • • , d, (2.8) 

where wi{t),--- ,Wdit) are independent standard Brownian motions on M. Thus 
it follows from (2.8) that the radial part of the d-dimensional Ornstein-Uhlenbeck 
process v(t) is given by 



i?'(t) : = V{MtW + --- + {vd{t)) 



2 



ae-'^'R'^ 



-,2at 



l\ (2.9) 



2a 



where R'^ is a Bessel process of dimension d starting from 0. Thus this motivates 
the following definition. 

Definition 2.2. Let (5 > 0, and let a and a be as above. If {R^{t))t>o is a Bessel 
process of dimension 5 starting from 0, then the process {R^{t))t>o defined by 

/ 2at _ 1 \ 

i?^(t) = ae-"*i?M— j (2.10) 

is called the generalized Bessel process of dimension 5 starting from with param- 
eters (a, cr). 

From the above definition, a straightforward calculation using (2.2) and (2.3) 
shows that the pdf and the transition probability density function of R^{t) are 
given by 

p?(x)=pf(O,x) = 2-i+V(i)-^r(0 ^5-iexp(-^) l[o,oo)(^) (2.11) 



and 

J 



Pt{xo,x) 



Pi^) ' T^^^rr ^exp 



e '^^xq 



^-1 



(x^ + e-2''*a;^) 



2p{t) 



xe '^^xq 

Is_,[—^],Xo>0. 



(2.12) 
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Remark 2.3. In [Eie], the author used the term generalized Bessel process to de- 
note the radial part of a d-dimensional Ornstein-Uhlenbeck process. Here we are 
adapting this terminology to the more general situation when 5 is not necessarily 
an integer. 

lfpt{x) denote the pdf of the Ornstein-Uhlenbeck process above, it is well-known 
that [N] 

/•OO 

\im pt{xo,x) = p^{xo)ptXo,x)dxo=Poo{x), (2.13) 



f— )-oo 

where Poo{x) is the pdf of the normal distribution A^(0, p(oo)). 

Proposition 2.4. For the transition probability density function of the generalized 
Bessel process in (2. 11), (2. 12), 

limpf{xo,x)= pi^{xo)pf{xo,x)dxo=pl^{x), (2.14) 

where 

pI{x) = 2-i+V(oo)-i+ir (0 ' x^-^exp (-^) • (2.15) 

Proof. The transition probability density function qi (xq , x) of the square of the 
generalized Bessel process above is related to pI{xq,x) by 

ql{xQ,x) = -^pl{^,^). (2.16) 

zyx 

Hence the explicit expression for ql{xQ,x) can be computed from (2.11) and (2.12). 
Note that the function g^(x) = 2~^x~'^/'^p^^[^/x) (resp. q^{xQ^x)) is related to the 
one in Eqn.(28) of [W] (resp. Eqn.(30) of [W]) by scaling of the time variable, the 
spatial variables and an overall scaling of the function itself. Hence it follows from 
[W] that 

/•OO 

lim ^f(a;o,x) = / ql^{xo)qlxQ,x)dxQ = qi^{x), (2.17) 



t— )-oo 

where 



5\-^ s 



qlix) = 2- v(oo)-^r i^-j x^-^exp [-^^^) ■ (2-18) 

Since p^{x) = 2xq^{x'^), it is immediate from (2.17) that liuit^ooPt{xo,x) = 
p^{x). Moreover, on using (2.15) and (2.17), we obtain 

/"OO 

pLixo)pt{xo,x)dxo=2x qL{yo)qf{yo,x^)dyo 
Jo (2-19) 

= 2xqlix^)=pUx)- 
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So this completes the proof. D 

Remark 2.5. (a) Alternatively, one could give a direct proof of the limiting be- 
haviour of pf (xq, x) for xq > by invoking the asymptotics I,y{z) r^ {z/2Y /r{u + l) 
as z — )■ 0. On the other hand, the fact that p^{x) is the stationary density of the 
generalized Bessel process can also be established in a direct way by making use of 
the Weber Sonine formula [Wat] 

x^e"^V^(7x)(ix 

valid for Re a > 0, Re(;U + v) > —1. We will leave the details to the interested 
reader. 

(b) In [W] , the author applies the method of separation of variables to the Fokker- 
Planck equation and identifies the cofficients of this equation with those of the 
Pearson's equation. In this way, eigenvalue problems of Sturm-Liouville type with 
reflecting boundary conditions at the end points are obtained. It is interesting to 
point out that as a result of the connection which we mentioned above, we have 
the expansion 

Pt(xo,x) 

where B{a, b) is the Beta function and 

1 d'^ 

L"(x) = — x-°e-^^(x°e-^) (2.22) 

are the Laguerre polynomials. 

For our construction in the rest of the paper, we will pick the normalization 
(a, a) = (1/2, 1). In this case, p^{x) is just the pdf of the Chi distribution xs while 
Pooix) is that of standard normal. Hence the corresponding generalized Bessel 
process and Ornstein-Uhlenbeck process have the desired properties. (See (1.1) 
above.) 



L.-C. LI 



3. The beta-Hermite processes. 



We begin by defining the product Markov process which will be used in our con- 
struction in the present section and the next. For this purpose, let X = {Xt)t>o and 
Y = {Yt)t>o be two statistically independent time-homogeneous Markov processes 
with transition probability density functions given by px{t,XQ,x) and PYit,yo,y) 
respectively. Then the product process X f^Y is defined to be the process where 
{X0Y)t = {Xt,Yt) fort > 0. Futpx(s,Y{t,{xo,yo),{x,y)) = px{t,xo,x)pY{t,yo,y), 
then it can be verified that the product process is again a Markov process with 
Pxi^Y{t,{xo,yo),{x,y)) as its transition probability density function. Clearly, we 
can iterate this construction and so we can define the product Markov process for 
any given number of statistically independent time-homogeneous Markov processes. 

3.1 The beta-Hermite processes and their eigenvalue distribution. 

Definition 3.1.1. The /3-Hermite process (J/3(t))t>o is the stochastic process on 
n X n Jacobi matrices Ji3{t) given by 



Mt) 



( 7^^^^'^ 


^R^-'^nt) 


•••\ 


^R^-'^^it) 


^^^(*) 


^R'-'^'it) ••. 





^R^-^^^it) 


• 



(3.1) 



•/ 



where the processes on the diagonal and the subdiagonal are statistically indepen- 
dent of each other. Here, f7i(t), • • • , f/„(t) are Ornstein-Uhlenbeck proceeses start- 
ing from with parameters (1/2, 1) and R^^{t) is the generalized Bessel process of 
dimension jf3 starting from with parameters (1/2, 1), j = 1, • • • , n — 1. 

Proposition 3.1.2. The f3-Hermite process is a matrix-valued diffusion process 
starting from, with transition probability density function 



n-l 



P{t,J,J) = 2^/?^ 



llpt{VPd,,y^a^)l[p^^ '^\^b,, 



2/3 6,), (3.2) 



i=l 



J = l 



^r'' 



with respect to Lebesgue measure dadh = dai ■ ■ ■ dondbi ■ ■ ■ dbn-i on M" x 
where di = Ju, hi = Ji,i+i, a^ = Ju and hi = Ji,i+i. Consequently, the joint density 
of the independent entries of Jp (t) is given by 



P{t,0,J)=Cnpp{t)-^ 



/3 r 



n-1 



^^ne-.^-p 



fc=i 



/3 



Mt) 



trJ' 



(3.3) 
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where 

2t-l^t + |n(n-l) 

Cn/3 = : 777V- (3-^) 

Finally, ifWh{a, b) is the joint density of the independent entries of the f3-Hermite 
ensemble in (1.1), then 

lim P{t, J,J)= [ Whid, b)P{t, J, J) dMb = Whia, b). (3.5) 

Proof. It is easy to see that the transition probabihty density functions of — ==i?'^(t) 
and -j^v{t) are given by ^/2|3pl{^/2|3 XQ,^/2|3x) and \/^pt{y/P xo,\fp x) respec- 
tively. Therefore, (3.2) is a consequence of the product construction of independent 
Markov processes. As P[J^(0) = 0] = 1, the pdf in (3.3) now follows from (3.2), 
(2.11), and the explicit formula for pt{x). Finally, the validity of (3.5) is due to 
(2.13) and (2.14) and this completes the verification. 

D 

Our next goal is to calculate the joint probability density function of the eigen- 
values of Jfi{t). In preparation, observe that the pdf of -j=Ui{t) can be expressed 
in terms of the pdf Poo(a;) of A^(0, 1). Indeed, this pdf is given by 



P<=o\ \ —7-r X 



(3.6) 



Similarly, it follows from (2.11) that the pdf of -^=W^{t) takes the form 



2 -^ / fl^2- 



V^pf (V^x) = ^'^'^,^. x^-^-^exp (-^) l[o,oo)(x) 



(3.7) 



^%- J^^ 



where p'^{x) is the pdf of the Chi distribution Xjp- Hence we have 

(n \ n-^ n n-1 



i=i j=i 



' -wj " 



(3.8) 



p{tr-h \^/W)'VW), 
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Now recall that a generic Jacobi matrix J = J {a, b) has nonzero entries on the 
subdiagonal and the eigenvalues of J are distinct. We will order the eigenvalues 
such that Ai > A2 > • • • A„ and denote by /i(l) > 0, /2(1) > 0, • • • , /„(1) > the 
first components of the normalized eigenvectors corresponding to the distinct eigen- 
values. Put A = (Ai, • • • , A„), /(I) = (/i(l), • • • , /n(l)). It is well-known that the 
map (p : J{a, b) — > (A, /(I)) is a diffeomorphism from the set of generic Jacobi ma- 
trices with positive subdiagonal entries to C^ x S*"" , where C^ = {{xi, • • • , x„) G 
M" I Xi > ■■■ > Xn}, and S*"""^ consists of vectors q = {qi,--- ,qn) on the unit 
sphere S"""^ such that qi > for all i (see, for example, [Dl]). Moreover, it follows 
from [DE],[D2] that 

T-rra — 1 I 

dadb= J.t=\/ dXda, (3.9) 



where 



4fi(i)---4f„_i(i) 



is the element of surface area in S*" 

Theorem 3.1.3. Under the /3-Hermite process, the vector X{t) = (Ai(t),--- ,A„(t)) 
of ordered eigenvalues of Jj^it) and the vector /(l,t) = (/i(l,t),--- ,/„(!, t)) of 
first components of normalized eigenvectors are independent. If A (A) is the Van- 
dermonde determinant, then 

(a) the joint pdf of the unordered eigenvalues of J^ (t) is given by 

Pn,{t, A) = C„, p(t)- t-^(-i) I A(A)|^ exp [-^ E ^') 



Cnp p(t)- t-T«("-i) eM-m{t, A)), 



where 



and 



JL r(i + f 



C, = (2vr)-t/?t+f "(-1) n -4 ^, (3.12) 



=ir 1 + ^ 



1 '^ 

W{t,X) = —-S2x^ -Y.^og\X, - X,\. (3.13) 

"^P^^) i=i i<j 

(b) the joint density 0/ /i(l, t), • • • , /„(1, t) with respect to the measure da on S"~ 
is given by 

r (l^\ n 
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Proof. Let (p be the map above and consider J^(t) in the domain of (p. Let Ai(t) > 
• • • > A„(t) be the (ordered) eigenvalues of J^(t) and let /i(l, t) > 0, •••/„(!, t) > 
be the first components of the normalized eigenvectors corresponding to the eigen- 
values. Because of (3.8), the independence of A(t) and f{l,t) follows as in the 
calculation in [DE] where we have to use (3.9) and the relation [D1],[DE] 

A(A) ^ n(^^ - ^^o = nr\":; - (3.15) 

The calculations leading to the assertions in the remaining parts of the proposition 
are also similar and so we skip the details. D 

Now, another way to describe generic Jacobi matrices J{a, b) is by using the 
spectral measure 

n 

i=i 
Indeed, it is well-known that the map ^p : J{a, b) — > ;U is a bijection from the set 
of generic n x n Jacobi matrices to the set of probability measures on M supported 
at n points. Let 

/x(t) = (^1 (*),-•• ,^„(t)), fi,it) = f,il,tf,l<j<n. (3.17) 

To describe the probability distribution of the vector fi{t), recall that the Dirichlet 
distribution Dir„_i(Qi, • • • ; q„) with parameters ai, • • • , a„ > is the distribution 
which has a density with respect to Lebesgue measure on M""-*^ given by (see [Wi] 
for more details) 

r(a,)---r(«n) n^/ i^(2/ir--,y.-i) O-is) 

where 

y^ = l-yi y„_i, (3.19) 

and where S is the simplex 

S=l{yi,--- , yn-i)\ %• > for all j, J] %• < 1 I . (3.20) 

Note that for n = 2, Diri(ai;Q2) is the Beta distribution Beta(ai,a2) which is 
supported on [0,1]. The following result is a straightforward consequence of the 
joint density of /i(l, t), • • • , /„(1, t) with respect to da and the basic properties of 
the Dirichlet distribution [Wi]. 
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Corollary 3.1.4. Under the f3-Hermite process, the vector of weights n{t) of the 
spectral measure associated with Jj3{t) follows the distribution Z)ir„_i(^, • • • ; ^). 
Hence the marginals are Beta distributions: 

/i,(t) ~ Beta (^, i!l_lM^ , 1 < J < n. (3.21) 

Moreover, for each 1 < k < n, 

EM,(*)~B.«»(f,^^). (3.22) 

3.2 Time-dependent semicircle law. 

We introduce the following scaling of J pit): 

j(;\t) = ^, (3.23) 



and let a)^ (t) = {Ja^ {t))kk, ^k {t) = {J p \t))k,k+i- The goal of this subsection is 
to study the large n behaviour of the spectral measure process 

n 

and the empirical eigenvalue process 



1 " 

;^E^.-W (3-25) 



i=l 

for each t > 0, where A^" (t), • • • , \n {t) are the eigenvalues of jJ^{t). We begin 
with a lemma. 

Lemma 3.2.1. As n ^ oo, 



arW^O, 6r(t)^y4^ (3.26) 

for each t > 0. 
Proof. For any t > and any r = 1, 2, • • • , it follows from (3.6) that 

1/2 / ^fl 9^ 



1 ^pW^ " 



2 /•"" _^ 

2t: \n(3 J J _^ 



(3.27) 
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from which it is clear that E[{aj^ {t)Y] — )■ as n — )■ oo. Therefore, aj^ (t) — > 
and hence a)^ (t) — > 0. Similarly, we obtain from (3.7) that 



E[{b^^\t)y 



nl3 
W) 



(n-fc)/3/2 



(n-fc)/3 



n/3 



p {(n-k)£\ _/q 
71(3 



^^^_k)p-l+r ^^^ 



„{n-k)^+r-\^-x^ ^^ 



n(3x^ 



dx 



(3.28) 



{n-k)li+r 



Y ( {n-k)l3 



From the asymptotics of the Gamma function, we have 



(n-fc)/3+r 
2 



Y [ (n-fc)/3 



n/3^ ^ 

T 



(3.29) 



as n — 7- oo. Hence the assertion 6^(t) — > y 2^ follows from (3.28) and (3.29). D 

We now turn to the analysis of the spectral measure process (/i")t>o. First, from 
the relations 

(3.30) 



,^dfi^{x) = {er,{J^^\t)fer), A; = !,•••, 

it is clear that the moments f^x'^dfi'^ix) are polynomials in the entries of Jj^'^t). 
As the entries on the diagonal and subdiagonal of Ji" (t) are independent random 
variables, it follows from (3.26) that as n — )■ 00, 



ia^-\t),b'r\t),---,af{t),bf{t))^. 



0, 



Pit) 



,0, 



Pit) 



(3.31) 



for each fixed value of j and each t > 0. Hence by the continuous mapping theorem, 
(3.30) and (3.31), we obtain 

x^'dfi'lix) A (ei,(j(°°)n*))^ei) 



as n — )■ 00 



where J^°°> (t) is the Jacobi operator on £2 given by 



/ 



J(°°))(t) 



\ 



n . ./p(*) 
" V 2 



J^ 



■■■) 



(3.32) 



(3.33) 
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Now the orthogonal polynomials corresponding to J^°°' (t) are defined by 

sm{n9) 



P'ix) 



and it is easy to check that 



smf 



X = \/2p{t) cos 6 



PLix)Pnix)dfit{x) = Smn, 



where 






(3.34) 



(3.35) 



(3.36) 



(3.37) 



Combining (3.32) and (3.37), we obtain the first part of the following theorem. 

Theorem 3.2.2. (a) For each t > 0, the sequence (/u")„>i converges weakly, in 
probability, to the probability measure pt defined in (3.36). 

(b) For each t > 0, the sequence {vt)n>i converges weakly, in probability, to the 
same probability measure fif 

Proof. We have already proved (a). In order to establish (b), it suffices to show 
that 



Thus dfit is the spectral measure of J^°"' (t) so that 

(ei,(j(°°)(t))'=ei)= f x^dptix) 



dLp{f^7,^t 



as n — )■ oo, 



(3.38) 



where d^p is the Levy-Prohorov metric on the space of (Borel) probability measures 
on M. For this purpose, introduce the distribution functions F^^ , F^^ corresponding 
to /i" and f " respectively. Then from the definitions of the Levy-Prohorov metric 
and the Levy distance between distribution functions, we have 



dLp{f4,Vt) < sup \F^r^{x) - F^^{x)\ 





< max 




and so it suffices to show that 


max 




p 

— 7- as n — 7- 



(3.39) 



(3.40) 
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As Yli=i f^j i^) has the Beta distribution (see (3.22)), the rest of the proof of 
identical to that of Theorem 5.4 in [BNR]. For the sake of completeness, we give 
the main steps. First of all, from the density of the Beta distribution, we have the 
moments of ^. ;^ /x^" (t): 



E 



E^r(o 



vj=i 



fc/3 



n/3 






(3.41) 



from which we see that 



E 



j=i 



k 
n 



By using (3.41) and its special case in (3.42), we find 

' k{n — k) 



E 






O 



n^ 



(3.42) 



(3.43) 



so that 



E^ 

fc=i 



E'-rc)-! 



.j=i 



oil 

n 



Hence for each e > 0, we obtain 



(3.44) 



P 



max 






> e 



< 



E^ 

fc=i 



E^r(^) 



i=i 



n 



> e 



<^-^E^ 



o 



fc=l 

1 



i=i 



n 



(3.45) 



from which (3.40) follows. 
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4. The beta-Laguerre and the beta-Wishart processes. 

The matrix model of the /3-Laguerre ensembles parametrized by a > —1 [DE] is 
defined schematically by 

1.., ,_ 1 .... , 



/ -^X{a+n)P -^Xl3{n-1) 



'"\ 



L 



I3,a 



(4.1) 



-^X(a+n-l)l3 ■^X(n-2)/3 

V : ••. ••• •••/ 

where the entries on the diagonal and the subdiagonal of the n x n matrix -L^,a 
are independent. Using the pdf of the Chi distribution, the joint density of the 
independent entries of L^.^ reads 



n-l 



Wdx,y) = d^,llxr-'^'^'-'e-^^'ll 



where x,- 



a+n-i+l)B-l -S.x'^ 
i=l i=l 

{Li3,a)ii, Ui = {Li3,a)i,i+i and 



(n-i)/3-l _£„2 



d 



2-1 , , . nal3 I 13 2 



nP 



rn-1 



n;^irff)n;=ir^'^^i±^'i 



(4.2) 



(4.3) 






In this section, we will present results analogous to those in Section 3 for two related 
matrix models: the one defined in (4.1) and an associated one consisting of matrices 

4.1 The beta Laguerre (Wishart) processes and the eigenvalue distribu- 
tion. 

Definition 4.1.1. The /3-Laguerre process {Lp^a{t))t>Q parametrized by a > —1 
is the stochastic process on n x n bidiagonal matrices Lp^it) given by 



/ li?(a+™)^(t) -i=i?("-l)^(t) 



Lp,ait) 



VP 



V?' 







-i=i?(«+"-)^(t) -^Ri 



V? 



v^ 



ra-2)/3 



it) 



\ 



...\ 



■••/ 



(4.4) 



where the processes on the diagonal and subdiagonal are statistically independent 
of each other, and where R^{t) is the generalized Bessel process of dimension 6 
starting from with parameters (1/2,1). Let Jp^ait) = L^^ait)'^ Lp^ait)-, then the 
process {Ji3,a{t))t>o is called the associated /3- Wishart process. 

As in the analogous case in Proposition 3.1.2, the following is immediate from 
the product construction of independent Markov processes, P[L^^a(0) = 0] = 1, 
(2.11) and (2.14). 
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Proposition 4.1.2. The j3-Laguerre process is the matrix-valued diffusion process 
starting from with transition probability density function 

n n — 1 

p{t, L, L) = r-'^ n^t"^""^^'^^(v^ *- v^ ^^) n ^f'^^^^^r^ y- ^^ ^^^ 

i=l j=l 

(4.5) 
with respect to Lebesgue measure dxdy = dxi • • • dx^dyi • • • dyn-i on M" x M"~ , 
where Xi = L^, iji = L^ i_|_i, Xi = La and bi = Lj j_|_i. Thus the joint density of the 
independent entries of Lp^a{t) is given explicitly by 

n Q^2 n — 1 Q 2 

T-,fx n T\ J /,N_I1££_£2_ TT (a+n-i+l)B-l -^^ TT (n-i)B-l --^^ i a n\ 

P{t,0,L) =dnfsp{t) 2 2 |J_x^^ ^ ^ "^ e ^^^ [[yl "^ e ^^^ . (4.6) 

i=l i=l 

IfWi{x,y) is the density in (4-1), we have 

lim P{t,L,L)= f Wiix,y)Pit,L,L)dxdy = Wiix,y). (4.7) 

We next calculate the joint probability density function of the eigenvalues of 
Jl3,a(t) = Lp^a{i)'^ Li3^a{i)- To do that, we have to first compute the pushforward 
of the measure P{t,0,L)dxdy under the map L ^^ J = L^ L which sends n x n 
bidiagonal matrices to n x n Jacobi matrices. To this end, observe that 

Pit L)= (^Y ' rr „(«+"-i+i)/3 ( ^O TT n(n-i)P y^ ] 



Therefore, if we let a^ = Ju = {L'^L)ii, bi = Ji^i+i = {L'^L)i^i^i, then in the 
notations of section 3 for Jacobi matrices, the calculations leading to the results in 
the next proposition are similar to the corresponding one in [DE] where one has to 
use (3.9), (3.15) and the fact that the Jacobian of the map L i— t- J = L^L is given 

by [DE] 

n-l 

2'^x„ n xl (4.9) 

i = l 

Theorem 4.1.3. Under the (3 -Wishart process, the vector \{t) = (Ai(t), • • • , A„(t)) 
of ordered eigenvalues of Jjs^ait) and the vector /(l,t) = (/i(l,t),--- ,/„(!, t)) 
of first components of normalized eigenvectors are independent. If A (A) is the 
Vandermonde determinant, then 
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(a) the joint pdf of the unordered eigenvalues of Jp^ait) is given by 

n f Q ^ \ 

^„",(^,A) = C„^,p(t)-^-f">(A)|^n\"^"^'^"-P (-^^Y.>^) ' (4-10) 
where 



, fl\ 4""-r4" " r(i + f 



(h) the joint density of fi{l,t),- ■ ■ , /„(1, t) with respect to the measure da on S_^ 
is given by 

rff 



n-l 



2""' / /.Vx- IIgf"'- (4-12) 



r(f)) .=1 



In view of (4.12), it follows that if fi(t) = (/ii(t),--- , finit)) is the vector of 
weights in the spectral measure of the spectral measure of Jp^ait)-, then the distri- 
butions of jj.j{t) and J2j=il^ji^) ^^6 ^Iso given by (3.21) and (3.22) respectively. 

4.2 Time-dependent Marchenko-Pastur law and quarter-circle law. 

We introduce the following scaling of L^^a(t): 

4:iw = ^' (4-13) 



and let x^^\t) = (L£(t)),,, y^^\t) = (L£(t)),,fc+i. 

We begin with the analog of Lemma 3.2.1 for the /3-Laguerre process. 

Lemma 4.2.1. As n ^ oo, 

xPit)^./^), y(")(t)A7^, (4.14) 

for each t > 0. 
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Proof. For any t > and any r = 1, 2, • • • , it follows from the pdf of R (t) that 

(71-fc)g 

2 ( -^^\ ^ 



2(2£WV roo 



n/3 ^ f{n-k)f3 + r 



(4.15) 



r. 2 



(ra-fc)/3 



as n — )■ oo where we have used the asymptotics of the Gamma function. Hence 
Vk v) — ^ \/p(t) as n — ;• oo. Similarly, 



(g + n — fc + l)/3 



2 r^i^^ ^ oo 

^^ fc ^ ^^ J r A«+»-fc+i)/? \ Jo V 2p(t) y 



r r 



(a+n-fc + l)/3+r 



r. 



\ n(3 J p ^ (a+n-fc + l);3 



(4.16) 
as n — )• oo and so we also have x)^ (t) — > yj p{t). D 



Now we introduce 



4"iw = (4:iw)^4:iw (4-17) 



and let al''{t) = {Jp^a{^))kk, o^ v) = ('^fi'ai'^))k,k+i- Then from the relations be- 
tween the entries of J^lit) and -Ll"^(t), we can deduce the following from Propo- 
sition 4.2.1 when we invoke the continuous mapping theorem. 

Corollary 4.2.2. As n ^ oo, 

(4.18) 
for each t > 0. 
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Let Xi (t),--- ,Xn{t) be the eigenvalues of JTa{t). We consider the spectral 
measure process 

n 

,(») 



^" = EM"ni)V)w (4-19) 

and the empirical eigenvalue process 

n 

-— ;7E\<")m (4-20) 






associated with ji"j (t) . From the relations 



^kA,,n(^\ — („ ( ji^) u\\f^ 



x'=d^rW = (ei,(4:JW)'ei), k = !,■■■, (4.21) 

it follows that the moments J^x^dii^{x) are polynomials in the entries of Jp^i^)- 

Thus it follows from Corollary 4.2.2 and the continuous mapping theorem that 

x'^d/.rix) A(ei,(4^»(t))'=ei) (4.22) 

as n — 7- oo, where Ja^ (t) is the Jacobi operator on £2 given by 

'p{t) p{t) ••• 

T(^))^.^ _ I P(0 2p(t) p{t) O--- I . . 

•^/^.^ ^^^ ~ I p{t) 2p{t) p{t) • • • I • ^ ^ ^> 



The next thing to do is to compute the spectral measure dp^t of JgJ{t)- For this 
purpose, we will make use of the method of Grosjean [G]. First of all, we introduce 
the (time-dependent) polynomials {P^(x)}„>o satisfying the recursion relations 

P^+i(x) = {x- 2p{t))Pi{x) - p{tfP^_,ix), n > 1 (4.24) 

and the initial conditions 

Po*(x) = l, Pl{x) = x-p{t). (4.25) 

Also, introduce the functions of the second kind 

Qi{z)= f ^Mdp,(x) (4.26) 

./» z — X 
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for n > and for 2; G C From (4.24) and (4.25) above, we have 

Q\{t) = -l + {z-p{t))Ql{z), (4.27) 

and 

QU,{z) = {z- 2p{t))Qi{z) - p{tfQi_,{z), n > 1, (4.28) 

where Q\){z) is the Stieltjes transform of the spectral measure dpf By making use 
of these relations, we obtain the continued fraction expansion 

Oo(^) = ,-,(t)'-F.(.) - ("■'''> 

where 

F\z) = P^'^' . (4.30) 

^-Mt) '^'^ ,,,, 

z - 2p{t) P^'^ 



z - 2p{t) 
But from the expression for F*{z), it is clear that 



z-2p{t) -F^{zy 
Solving, we obtain 



F\z) = ^ ''^^^ ^^, , . (4.31) 



F\z) = ^-Mt)-^/{z-2p{t)f-Ap{tY ^ ^^_^2^ 

where the branch of the square root is the one which tends to the positive square 
root of (a; — 2p(t))^ — 4p(t)^ when z tends to x G (4p(t),oo). Therefore, when we 
substitute (4.32) into (4.29), the result is 

dixt{x) 1 



/ 



z-x L + y(^,-2p{t)f-Ap{tf 



(4.33) 



We next introduce the Fourier transform of dpt- 

pt{x)= [ e'^^'dfitis). (4.34) 

Then from (4.33), its Laplace transform is given by 



C{pt){p)= / e-^PMx)dx 
Jo 

djitis) 



-ip — s 
i 



(4.35) 



f-iV(^P + 2p(t))^-4/>(t)2 
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for Re (p) > 0. Thus by the inversion theorem for Laplace transform, we have 

' 0, X > 

i/2t(0+), x = 



1 

2^ 



^xp 



tp 
(7 — ioo "2" 



i./ii^T2i^w^wy' 



dp 



(4.36) 



_ fit{x), X > 0, 

where o" is a real number which is greater than the real parts of the singularities of 
C{jlt){p) in the j»-plane. Now make the change of variable p = iu in the integral on 
the left hand side of (4.36), this gives 

f 0, X > 



2^ 



^ + yiu-2p{t))^-Apity 



i/2t(0+), x = 
_ Jlt{x), X > 0, 



(4.37) 



where the path of integration is now a horizontal line below the singularities of 
the integrand in the n-plane. But from the above expression, it is easy to show 
that the denominator of the integrand vanishes only at n = 0. Thus the set of 
singularities of the integrand coincides with the branch cut [0, 4p(t)] of the function 
\/{u — 2p{t))'^ — 4/9(t)2 and hence we can take a = — e, where e > is a small 
number. But the fact that dpt is a real measure means that /2((— x) = /if(x) for all 
X. Consequently, 

ds, (4.38) 



Ux) = - 



^ + \^u{u-Ap(t)) 



u=s — ze 



for all X G M where /it(0) = ^(/lt(0+) + /lt(0— )). But from the definition of 
\/u{u — 4/9(t)), we have 



lim ^/u{u^^4p(fy) \u=s±ie-- 

e->0+ 



f -^s{s-4p{t)), sG(-oo,0] 
±i^si4pit) - s), sG(0,4p(t)) 



(4.39) 



I ^ys{s-4p{t)), s G [4p(t),oo). 

Therefore, when we take the limit as e — )■ 0+ in (4.38), the only contribution to the 
integral comes from [0, 4p(t)] and we find 

.4p(t) -I 



Mx) 







2Trp{t) 



4p{t) 



ds. 



Hence we can now conclude that 



dpt{x) 



1 



27rp(t) 
and so we have the following result. 



Ap{t) - X 



l(o,4p(f)) dx 



(4.40) 



(4.41) 
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Theorem 4.2.3. (a) For each t > 0, the sequence {fi^)n>i converges weakly, in 
probability, to the probability measure fit defined in (4-4i)- 

(b) For each t > 0, the sequence (^'")n>i converges weakly, in probability, to the 
same probability measure fif 

We now return to the /3-Laguerre process itself. Note that although we have 
Theorem 4.2.3 available to us, however, it is not hard to see that it is not possible 
to deduce the corresponding result for the /3-Laguerre process. Our study of the /3- 
Laguerre process will be based on the following transformation which is well-known 
in numerical linear algebra [GK]. Suppose B is the bidiagonal matrix 



B 



/Xi 



Vo 



yi 

X2 



O \ 



y2 



Vn-l 



(4.42) 



with singular value decomposition B = UT,V'^ , where U = (ui, • • • , n„) and V = 
(fi, • • • ,Vn) are orthogonal, and S = diag (o"i, • • • , cr„) is the diagonal matrix whose 
diagonal entries are the singular values of B. Then the eigenvalues of the 2n x 2n 
symmetric tridiagonal matrix 

yi 





Xi 



S 



Xi 





yi 



Jbr, 



(4.43) 



V O a;„ / 

, iy„, — 0"i, • • • , — (T„ and the first components of the corresponding normal- 

l)/^/2, v^{l)/V2, 



are ci , • • • ,a. 

ized eigenvectors are given by t;i(l)/\/2, • • • ,f„(l)/-v/2, vi{\)/\/2, ••• , f„(l)/\/2 

respectively. In view of this, it suffices to study the process (5'^,a(i))t>0) where 



Sp,a{i) 



( 



V ; 



l^(a+n)^(i) 



V^ 







•••\ 



-^i?("-i)^(t) 



V? 



^i?("-l)^(t) 



(4.44) 



•••/ 



is obtained from (4.4) by applying the above transformation. (The full justification 
is in Theorem 4.2.4 below.) We introduce the following scaling of Sp^ai^)'- 



s',:iit) 



n 



(4.45) 
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(n), 



»l 



corresponding to (4.13). By our discussion above, if a\ (t), ■ ■ ■ ,an (t) denote the 
singular values of L^^a(t), then the eigenvalues of Sg'l{t) are given by cr}" (t), • • • , 
cjn (t), — cJi (t), • • • , —an{t). Now wc iutroduce the spectral measure process 



-2n 



Mt 



E^i(*)fc")m+'^-.(") 



(t) 



j=i 



and the empirical eigenvalue process 



■-271 



1 " 

7 I 5 (n)/,s + 5 (n) 



i=l 



'w 



associated with (5'i"^(t))t>05 where 



(4.46) 



(4.47) 



/I"(t) = -/j(l,t)2, j = l,...,n. 



(4.48) 



By using Lemma 4.2.1, and following the same procedure in Section 3.2, we 
obtain 

x^d-fil'^ix) A (ei,(^/2j(°°)) (*))'= ei) = [ x''djit{x), (4.49) 

where ji°°' is the Jacobi operator in (3.33) and 



VMt)-x ^ 

dnt{x) = z — 1 



27rp(t) [-V'^W'^VpW 1 



[x] dx. 



Thus we have proved 



(4.50) 



dLp{f^t"',IJ't 



as n — )■ oo. 



(4.51) 



We next show that 



dLp{lJif^, I't"') — > as n -^ oo. 



(4.52) 



Here the idea is also the same as before. First we calculate the distribution of the 
vector of weights 

m = i'fiiit),--- ,rnit)) (4.53) 

which is a generalized Dirichlet distribution (we can think of the one defined in 
(3.18) as the standard one) supported on the simplex 



n-l 



3(1/2) = {{y,,--- ,yn-i) 



y, > for all j, ^ y, < 1/2 



> . 



(4.54) 
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Then for each 1 < k < n, we find 

k 



^M^(t)~Beta(°'^' 



i=i 



fc/3 {n-k)(3 



(4.55) 



where Beta^ ''^' denotes the generahzed Beta distribution supported on [0, ^]. Since 
we have 

k 



d{jir,vn<^^^ 



E/^^W-^' (4-56) 

i=i 

the analysis proceeds as in Section 3.2. Consequently, when we combine (4.51) and 

(4.52), we conclude that 



di^r,jit 



as n 



oo. 



(4.57) 



We are now ready to state the main result for the /3-Laguerre process. To this 
end, introduce the empirical singular value process 



1 " 

— > 5 (n) 



(<) 



(4.58) 



("), 



associated with {L\ ^(i))t>o- Also, let 



y/4p{t) - X^ ^ 

^^^("^ = npit) ^ro 



[0,2v/pW ] 



(x) dx. 



(4.59) 



We will consider ly^ and pt as measures on R_(_ = [0, cx)). 

Theorem 4.2.4. For each t > 0, the sequence (^'")n>i converges weakly, in prob- 
ability, to the probability measure fit- 

Proof. We will deduce this result from (4.57). For this purpose, it is more conve- 
nient to use the bounded Lipschitz metric, which is equivalent to the Levy-Prohorov 
metric [Dud]. In order to write down the expression for this metric, denote by 
BL{S) the class of bounded functions / : S — > M on a complete metric space S 
which are Lipschitz. For / G BL{S), define the Lipschitz semi-norm 

\f(.x)-f{y)\ 



L(S) = SUp^^y 



(4.60) 



d{x,y) 

and put \\f\\BL(s) = II/I|l(s) + II/I|l<»(s) where ||/||l<»(5) is the sup-norm. Then 
||/||bl(5) is a norm and {BL{S), \\-\\bl(s)) is a Banach space. With these notations, 
the bounded Lipschitz distance between the two measure z/|" and Jit is given by 



dBL{R){j^t"',l^t) =SUPf 



eBLi 



fdv^ 



fdpt 



(4.61) 
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where the supremum is taken over 



BLi{R) = {f€BL{R)\ ||/bL(R) < !}• 
Now let BLl{R) = {/ G BLi{R) \ /is even}, then clearly 



(4.62) 



dBL{R)ii't",l^t)> sup f^smvL) fdi^t" - fdl^t ■ (4.63) 

Jr Jr 

But for / G BLl (M) , it follows from the definition of the two measures that 

/ fdu^- =^J2f{af^(t)) = I fdu^, (4.64) 

while 



j=i 



/ fdjit = I 

JR JR- 



f{x)dfif 



(4.65) 



As we can identify the space BLf{R) with BLi{R^), when we combine (4.63)-(4.65), 
the result is 



dBL(R)('^t",^t) >SUp/gs^^(jj^) / fdul 

Jr+ 

= dBL(R+){j^t->P't)- 



fdnt 



(4.66) 



Hence it follows from (4.57) and (4.66) that d_BL(R+ )('^r5 Mt) — ^ as n — )■ oo. 



Remark ^.2.5. If we let 



D 






(4.67) 



then following the same argument as in the proof of Theorem 4.2.4, we also have 

dBL{R+){fit,fJ't) — ^0 asn-^oo. (4.68) 

The measure valued process (/Z")t>o, however, is the spectral measure process of 
the square root ( d J^lit) I of the scaled /3-Wishart process, where Jg^i^) ^^ 



t>o 



defined in (4.17). 
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